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Abstract 



It is a deceptively simple question to ask how acoustic disturbances 
propagate in a non-homogeneous flowing fluid. Subject to suitable re- 
strictions, this question can be answered by invoking the language of 
Lorentzian differential geometry. I begin this paper with a pedagogical 
derivation of the following result: If the fluid is barotropic and inviscid, 
and the flow is irrotational (though possibly time dependent), then the 
equation of motion for the velocity potential describing a sound wave is 
identical to that for a minimally coupled massless scalar field propagating 
in a (3 + l)-dimensional Lorentzian geometry 



The acoustic metric g^(t, x) governing the propagation of sound depends 
algebraically on the density, flow velocity, and local speed of sound. Even 
though the underlying fluid dynamics is Newtonian, non-relativistic, and 
takes place in flat space plus time, the fluctuations (sound waves) are gov- 
erned by an effective (3 + l)-dimensional Lorentzian space-time geometry. 
This rather simple physical system exhibits a remarkable connection be- 
tween classical Newtonian physics and the differential geometry of curved 
(3 + 1) -dimensional Lorentzian spacetimes, and is the basis underlying a 
deep and fruitful analogy between the black holes of Einstein gravity and 
supersonic fluid flows. Many results and definitions can be carried over di- 
rectly from one system to another. For example, I will show how to define 
the ergosphere, trapped regions, acoustic apparent horizon, and acoustic 
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event horizon for a supersonic fluid flow, and will exhibit the close rela- 
tionship between the acoustic metric for the fluid flow surrounding a point 
sink and the Painleve-Gullstrand form of the Schwarzschild metric for a 
black hole. This analysis can be used either to provide a concrete non- 
relativistic analogy for black hole physics, or to provide a framework for 
attacking acoustics problems with the full power of Lorentzian differential 
geometry. 

1 Introduction 

In 1981 Unruh developed a way of mapping certain aspects of black hole physics 
into problems in the theory of supersonic acoustic flows jD. The connection 
between these two seemingly disparate systems is both surprising and powerful, 
and has been independently rediscovered several times over the ensuing decade 
and a half [^J. Over the last six years, a respectable body of work has been 
developed using this analogy to investigate micro-physical models that might 
underly the Hawking radiation process from black holes (or acoustic holes - 
"dumb holes"), and to investigate the extent to which the Hawking radiation 
process may be independent of the physics of extremely high-energy trans- 
Planckian modes §,g§§0§|0g|l![g,@. 

In this paper, I wish to take another look at the derivation of the relationship 
between curved spacetimes and acoustics in flowing fluids, to provide a pedagog- 
ically clear and precise derivation using a minimum of technical assumptions, 
and to develop the analogy somewhat further in directions not previously envis- 
aged. In particular, I will show how to define the notions of ergo-region, trapped 
regions, acoustic apparent horizons, and acoustic event horizons (both past and 
future) for supersonic fluid flows, and show that in general it is necessary to 
keep these notions distinct. 

As a particular example of a simple model exhibiting such behaviour I write 
down the acoustic metric appropriate to a draining bathtub [(2 + 1) dimensions], 
and the equivalent vortex filament sink [(3 + 1) dimensions]. 

I shall further show that the relationship between the Schwarzschild geome- 
try and the acoustic metric is clearest when the Schwarzschild metric is written 
in the Painleve-Gullstrand|] form ]l|, |l|, |l7|, p], [l9|, and that while the re- 
lationship is very close it is not exact. (It is in fact impossible to obtain an 
acoustic metric that is identical to the Schwarzschild metric, the best that one 
can achieve is to obtain an acoustic metric that is conformally related to the 
Schwarzschild metric.) If all one is interested in is either the Hawking temper- 
ature or the behaviour in the immediate region of the event horizon then the 
analogy is much closer in that the conformal factor can be neglected. 

For an arbitrary steady flow the "surface gravity" (mutatis mutandis, the the 
Hawking temperature) of an acoustic horizon will be shown to be proportional 

1 This is also often called the Lemaitre form of the Schwarzschild metric. 
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to a combination of the normal derivative of the local speed of sound and the 
normal derivative of the normal component of the fluid velocity at the horizon. 
In general, the "surface gravity" is 

1 d{c 2 -vl) d(c-v ± ) 
9H= 2 dn =C ^^- W 

(This generalizes the result of Unruh |j to the case where the speed of sound 
is position dependent and/or the acoustic horizon is not the null surface of the 
time translation Killing vector. This result is also compatible with that deduced 
for the solid-state black holes of Reznik [Q , and with the "dirty black holes" 
of©.) 

Finally, I shall show how to formulate the notion of a static (as opposed 
to merely stationary) acoustic metric and exhibit the constraint that must be 
satisfied in order to put the acoustic metric into Schwarzschild coordinates. I 
point out that while this is a perfectly acceptable and correct mathematical 
step, and a perfectly reasonable thing to do in general relativity, it is (I claim) 
a good way to get confused when doing acoustics — from the Newtonian view 
underlying the equations of fluid motion the Schwarzschild coordinate system 
corresponds to a very peculiar way of synchronizing (or rather, de-synchronizing) 
your clocks. 

To begin the discussion: I address the deceptively simple question of how 
acoustic disturbances propagate in a non-homogeneous flowing fluid. It is well 
known that for a static homogeneous inviscid fluid the propagation of sound 
waves is governed by the simple equation j2l|, 22 , 23 [m| 

dU = c 2 \7 2 yj. (2) 

(Here c = speed of sound.) Generalizing this result to a fluid that is non- 
homogeneous, or to a fluid that is in motion, possibly even in non-steady motion, 
is more subtle than it at first would appear. 

An important aspect of this paper is to provide a pedagogical proof of the 
following theorem: 



Theorem If a fluid is barotropic and inviscid, and the flow is irrotational 
(though possibly time dependent) then the equation of motion for the velocity 
potential describing an acoustic disturbance is identical to the d'Alembertian 
equation of motion for a minimally coupled massless scalar field propagating in 
a (3 + 1) -dimensional Lorentzian geometry 

A-0 = -Ld^ 9^ dui/i) = 0. (3) 

Under these conditions, the propagation of sound is governed by an acous- 
tic metric — g^,v{t,x). This acoustic metric describes a (3 + 1) -dimensional 
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Lorentzian (pseudo-Riemannian) geometry. The metric depends algebraically 
on the density, velocity of flow, and local speed of sound in the fluid. Specifically 



c 



~(c 2 -v 2 ) 



(4) 



(Here / is the 3x3 identity matrix.) In general, when the fluid is non- 
homogeneous and flowing, the acoustic Riemann tensor associated with this 
Lorentzian metric will be nonzero. 

It is quite remarkable that even though the underlying fluid dynamics is 
Newtonian, nonrelativistic, and takes place in flat space plus time, the fluctu- 
ations (sound waves) are governed by a curved (3 + l)-dimensional Lorentzian 
(pseudo-Riemannian) spacetime geometry. 

For practitioners of general relativity this paper describes a very simple and 
concrete physical model for certain classes of Lorentzian spacetimes, including 
black holes. On the other hand, the discussion of this paper is also potentially 
of interest to practitioners of continuum mechanics and fluid dynamics in that 
it provides a simple concrete introduction to Lorentzian differential geometric 
techniques. 

2 Fluid Dynamics 

2.1 Fundamental equations 

The fundamental equations of fluid dynamics |3[ |24| are the equation 

of continuity 

dtp + V • (p v) = 0, (5) 

and Euler's equation 

UD — * 

p =p[d t v + (v-V)v\=F. (6) 
at 

I start the analysis by assuming the fluid to be inviscid (zero viscosity) , with the 
only forces present being those due to pressure, plus Newtonian gravity, and with 
the inclusion of any arbitrary gradient-derived and possibly even time-dependent 
externally-imposed body force. Then 

F = -Vp - pV(j> - pV$. (7) 

Here (j> denotes the Newtonian gravitational potential, while $ denotes the po- 
tential of the external driving force (which may in fact be zero) j^] 

2 These two terms are lumped together without comment in fa|{, and are neglected in In. 
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Via standard manipulations the Euler equation can be rewritten as 

d t v = vx (V xv) - -Wp - V Qw 2 + <f) + . (8) 

Now take the flow to be vorticity free, that is, locally irrotational^ Introduce 
the velocity potential ip such that v — — V^, at least locally.^] If one further 
takes the fluid to be fraroiropjcQ (this means that p is a function of p only), it 
becomes possible to define 

f p dp' 1 
hip) = I -t-tt; so that Vh = - Vp. (9) 
Jo Pi?) P 

Thus the specific enthalpy, hip), is a function of p only. Euler 's equation now 
reduces to 

-d t ip + h + i(W>) 2 + + $ = 0. (10) 

This is a version of Bernoulli's equation in the presence of external driving 
forces. 



2.2 Fluctuations 

Now linearize these equations of motion around some assumed background 
[p ,Po,4>o)- Set p = p + epi + 0{e 2 ), p = p + epi + 0(e 2 ), and ip = 
t/jo + eipi + 0(e 2 ). The gravitational potential 4>, and driving potential i>, 
are taken to be fixed and external^. Sound is defined to be these linearized 
fluctuations in the dynamical quantities. Please note that this is the standard 
definition of sound and more generally of acoustical disturbances. In principle, 
of course, one is really interested in solving the complete equations of motion 
for the fluid variables ip,p,ip). In practice, it is both traditional and extremely 

3 The irrotational condition is automatically satisfied for the superfluid component of phys- 
ical superfluids. This point has been emphasised by Comer [E5|, who has also pointed out that 
in superfluids there will be multiple acoustic metrics (and multiple acoustic horizons) corre- 
sponding to first and second sound. Even for normal fluids, vorticity free flows are common, 
especially in situations of high symmetry. 

4 It is sufficient that the flow be vorticity free, V X v = 0, so that velocity potentials exist 
on an atlas of open patches — this enables us to handle vortex filaments, where the vorticity 
is concentrated into a thin vortex core, provided we do not attempt to probe the vortex core 
itself. I do not need tji to be globally defined. 

5 An unstated assumption of this type is implicit, though not explicit, in the analysis 
of reference [h|. On the other hand, reference B explicitly makes the stronger assumption 
that the fluia is isentropic. (That is, the specific entropy density is taken to be constant 
throughout the fluid.) This is a stronger assumption than is actually required, and the weaker 
barotropic assumption used here is sufficient. In particular, the present derivation also applies 
to isothermal perturbations of an isothermal fluid. 

e Fixed means that I do not allow back-reaction to modify the gravitational or driving 
potentials. Fixed does not necessarily mean time independent, as I explicitly wish to allow 
the possibility of time dependent external driving forces. 



5 



useful to separate the exact motion, described by the exact variables, (p,p,ip), 
into some average bulk motion, (po 7 po,ipo), plus low amplitude acoustic distur- 
bances, (epi,epi,eipi). See, for example |2f], |2^, |2^, [2^ |. 

Since this is a subtle issue that I have seen cause considerable confusion in 
the past, let me be even more explicit by asking the rhetorical question: "How 
can we tell the difference between a wind gust and a sound wave?" The answer 
is that the difference is to some extent a matter of convention — sufficiently low- 
frequency long- wavelength disturbances (wind gusts) are conventionally lumped 
in with the average bulk motion. Higher-frequency, shorter-wavelength distur- 
bances are conventionally described as acoustic disturbances. If you wish to be 
hyper-technical, we can introduce a high-pass filter function to define the bulk 
motion by suitably averaging the exact fluid motion. There are no deep physical 
principles at stake here — merely an issue of convention. 

The place where we are making a specific physical assumption that restricts 
the validity of our analysis is in the requirement that the amplitude of the 
high-frequency short-wavelength disturbances be small. This is the assump- 
tion underlying the linearization programme, and this is why sufficiently high- 
amplitude sound waves must be treated by direct solution of the full equations 
of fluid dynamics. 

Linearizing the continuity equation results in the pair of equations 

dtPo + V • (po v ) = 0, (11) 
d t pi + V • (pi v + po vi) = 0. (12) 

Now, the barotropic condition implies 

h{p) = h( Po + m + 0(e 2 )) = + e — + 0(e 2 ). (13) 

Po 

Use this result in linearizing the Euler equation. We obtain the pair 

-dtik + ho + \ (V^ ) 2 + + $ = 0. (f 4) 

-fit^i + — - v • VtAi = 0. (15) 
Po 

This last equation may be rearranged to yield 

Pi = Po{d t ^\ + v a ■ Vipi). (16) 
Use the barotropic assumption to relate 

dp dp 

Pl = dp Pl = dp Po + ^° ' ( 17 ) 

Now substitute this consequence of the linearized Euler equation into the lin- 
earized equation of continuity. We finally obtain, up to an overall sign, the wave 



G 



equation: 

(dp 
q- Po (d t tpi + v ■ Vtpi) 

+ V • (^p Vipi -^Povo (dtipi + v • V^i)j = 0. (18) 

This wave equation describes the propagation of the linearized scalar potential 
ipi. Once ipi is determined, equation ( fl6| ) determines pi, and equation (|l^) then 
determines p\ . Thus this wave equation completely determines the propagation 
of acoustic disturbances. The background fields p$, po and vq = — VV'o, which 
appear as time-dependent and position-dependent coefficients in this wave equa- 
tion, are constrained to solve the equations of fluid motion for an externally- 
driven, barotropic, inviscid, and irrotational flow. Apart from these constraints, 
they are otherwise permitted to have arbitrary temporal and spatial dependen- 



Now, written in this form, the physical import of this wave equation is 
somewhat less than pellucid. To simplify things algebraically, observe that the 
local speed of sound is defined by 



c- 2 = 



dp 
dp 

Now construct the symmetric 4x4 matrix 



Po 



(19) 



(20) 



(Greek indices run from 0-3, while Roman indices run from 1-3.) Then, intro- 
ducing (3 + 1) -dimensional space-time coordinates — x^ = (t; x l ) — the above 
wave equation ([l8]) is easily rewritten as 

W*" 9^0 = 0. (21) 

This remarkably compact formulation is completely equivalent to equation ( fL8| ) 
and is a much more promising stepping-stone for further manipulations. The 
remaining steps are a straightforward application of the techniques of curved 
space (3 + l)-dimensional Lorentzian geometry. 



3 Lorentzian Geometry 



In any Lorentzian (that is, pseudo-Riemannian) manifold the curved space 
scalar d'Alcmbertian is given in terms of the metric g^vit, x) by (see, for exam- 
ple, |2§ [27j [2§, ]29[ |30|) 

AiP = -^=0^ {^—g gT M . (22) 



-fJ 
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The inverse metric, g^itjx), is pointwise the matrix inverse of g/j, u (t,x), while 
g = det(<7 Ml /). Thus one can rewrite the physically derived wave equation 
in terms of the d'Alembertian provided one identifies 



-9 9 



f 



This implies, on the one hand 

det(/H = {^—gf 



= 9- 



(23) 



(24) 



On the other hand, from the explicit expression (p0|), expanding the determinant 
in minors 



P% 



Thus 



9 = A 
c z 



5-4 

c 



We can therefore pick off the coefficients of the inverse acoustic metric^] 

I 1 



f(t,Z) = — 
Poc 



(c 2 <5^ - v l vl) 



(25) 
(26) 

(27) 



We could now determine the metric itself simply by inverting this 4x4 matrix. 
On the other hand, it is even easier to recognize that one has in front of one an 
example of the Arnowitt-Deser-Misner split of a (3 + l)-dimensional Lorentzian 
spacetime metric into space + time, more commonly used in discussing initial 
value data in Einstein's theory of gravity — general relativity. (See, for example, 
[p8| pp 505-508.) The acoustic metric is 



9p,v 



Po 

c 



-I'n 



iv . u %] 

Equivalently, the acoustic interval can be expressed as 



(28) 



ds 2 



g^ v dx^ dx v 



Pa 



c 2 dt 2 + (dx l - v l dt) 5 i:j {dx 3 - v 3 dt) 



A few brief comments should be made before proceeding: 

• Observe that the signature of this metric is indeed (- 
should be to be regarded as Lorentzian. 



(29) 



as it 



7 There is a minor typo, a missing factor of c, in [E 
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• It should be emphasized that there are two distinct metrics relevant to 
the current discussion: 

— The physical spacetime metric is just the usual flat metric of Minkowski 
space 

r7 MV = (diag[-cf igW , 1, 1, (30) 

(Here cu g ht = speed of light.) The fluid particles couple only to the 
physical metric rj^. In fact the fluid motion is completely non- 
relativistic — ||«o|| < ciight- 

— Sound waves on the other hand, do not "see" the physical metric at 
all. Acoustic perturbations couple only to the acoustic metric g^ u . 

The geometry determined by the acoustic metric does however inherit 
some key properties from the existence of the underlying flat physical 
metric. 

• For instance, the topology of the manifold does not depend on the par- 
ticular metric considered. The acoustic geometry inherits the underlying 
topology of the physical metric — 3? 4 — with possibly a few regions ex- 
cised (due to imposed boundary conditions). 

• Furthermore, the acoustic geometry automatically inherits the property 
of "stable causality" g Note that 

gT (V„t) (V„t) = < 0. (31) 

This precludes some of the more entertaining causality-related pathologies 
that sometimes arise in general relativity. (For a discussion of causal 
pathologies, see for example [pl[ ). 

• Other concepts that translate immediately are those of "ergo-region" , 
"trapped surface" , "apparent horizon" , and "event horizon" . These no- 
tions will be developed fully in the following section. 

• The properly normalized four-velocity of the fluid is 

This is related to the gradient of the natural time parameter by 

V^= (1,0,0,0); V"t = ~^J=-- —. (33) 

p c Jp c 
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Thus the integral curves of the fluid velocity field are orthogonal (in the 
Lorentzian metric) to the constant time surfaces. The acoustic proper 
time along the fluid flow lines (streamlines) is 



and the integral curves are geodesies of the acoustic metric if and only if 
Poc is position independent. 

• Observe that in a completely general (3 + l)-dimensional Lorentzian ge- 
ometry the metric has 6 degrees of freedom per point in spacetime. (4x4 
symmetric matrix 10 independent components; then subtract 4 coordi- 
nate conditions). In contrast, the acoustic metric is more constrained. Be- 
ing specified completely by the three scalars ipo{t,x), po(t,x), and c(t,x), 
the acoustic metric has at most 3 degrees of freedom per point in space- 
time. The equation of continuity actually reduces this to 2 degrees of 
freedom which can be taken to be an d c(t, x). 

• A point of notation: Where the general relativist uses the word "station- 
ary" the fluid dynamicist uses the phrase "steady flow" . The general- 
relativistic word "static" translates to a rather messy constraint on the 
fluid flow (to be discussed more fully below). 

• Finally, I should add that in Einstein gravity the spacetime metric is 
related to the distribution of matter by the non-linear Einstein-Hilbert 
differential equations. In contrast, in the present context, the acoustic 
metric is related to the distribution of matter in a simple algebraic fashion. 

4 Ergo-regions, trapped surfaces, and acoustic 
horizons 

Let's start with the notion of an ergo-region: Consider integral curves of the 
vector = (d/dt)" = (1,0,0,0)^. (If the flow is steady then this is the 
time translation Killing vector. Even if the flow is not steady the background 
Minkowski metric provides us with a natural definition of "at rest".) ThenF] 



This changes sign when \\v\\ > c. Thus any region of supersonic flow is an 
ergo-region. (And the boundary of the ergo-region may be deemed to be the er- 
gosphere.) The analogue of this behaviour in general relativity is the ergosphere 

^Henceforth, in the interests of notational simplicity, I shall drop the explicit subscript 
on background field quantities unless there is risk of confusion. 




(34) 



9fiU (d/atr (d/dty = 9tt 




(35) 
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surrounding any spinning black hole — it is a region where space "moves" with 
superluminal velocity relative to the fixed stars [^8[ ^9], pPf . 

A trapped surface in acoustics is defined as follows: take any closed two- 
surface. If the fluid velocity is everywhere inward-pointing and the normal 
component of the fluid velocity is everywhere greater than the local speed of 
sound, then no matter what direction a sound wave propagates, it will be swept 
inward by the fluid flow and be trapped inside the surface. The surface is 
then said to be outer-trapped. (For comparison with the usual situation in 
general relativity see j2^, pages 319-323] or pages 310-311].) Inner-trapped 
surfaces (anti-trapped surfaces) can be defined by demanding that the fluid flow 
is everywhere outward-pointing with supersonic normal component. It is only 
because of the fact that the background Minkowski metric provides a natural 
definition of "at rest" that we can adopt such a simple definition. In ordinary 
general relativity we need to develop additional machinery, such as the notion 
of the "expansion" of bundles of ingoing and outgoing null geodesies, before 
defining trapped surfaces — that the above definition is equivalent to the usual 
one follows from the discussion on pages 262-263 of Hawking and Ellis [^9| . The 
acoustic trapped region is now defined as the region containing outer trapped 
surfaces, and the acoustic (future) apparent horizon as the boundary of the 
trapped region. (We can also define anti-trapped regions and past apparent 
horizons but these notions are of limited utility in general relativity.) 

The event horizon (absolute horizon) is defined, as in general relativity, by 
demanding that it be the boundary of the region from which null geodesies 
(phonons) cannot escape. This is actually the future event horizon. A past 
event horizon can be defined in terms of the boundary of the region that cannot 
be reached by incoming phonons — strictly speaking this requires us to define 
notions of past and future null infinities, but I will simply take all relevant 
incantations as understood. In particular the event horizon is a null surface, 
the generators of which are null geodesies. 

In all stationary geometries the apparent and event horizons coincide, and 
the distinction is immaterial. In time-dependent geometries the distinction is 
often important. When computing the surface gravity I shall restrict attention 
to stationary geometries (steady flow) . In fluid flows of high symmetry, (spher- 
ical symmetry, plane symmetry) the ergosphere may coincide with the acoustic 
apparent horizon, or even the acoustic event horizon. This is the analogue of 
the result in general relativity that for static (as opposed to stationary) black 
holes the ergosphere and event horizon coincide. 

5 Vortex geometries 

As an example of a fluid flow where the distinction between ergosphere and 
acoustic event horizon is critical consider the "draining bathtub" fluid flow. I 
model a draining bathtub by a (2 + 1) dimensional flow with a sink at the origin. 
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The equation of continuity implies that for the radial component of the fluid 
velocity we must have 

p v r cx -. (36) 



r 



In the tangential direction, the requirement that the flow be vorticity free (apart 
from a possible delta-function contribution at the vortex core) implies, via 
Stokes' theorem, that 



«*«-. (37) 
r 

On the other hand, assuming conservation of angular momentum (this places 
a constraint on the external body forces by assuming the absence of external 
torques) implies the slightly different constraint 

pv { oc-. (38) 
r 

Combining these constraints, the background density p must be constant (position- 
independent) throughout the flow (which automatically implies that the back- 
ground pressure p and speed of sound c are also constant throughout the fluid 
flow). Furthermore for the background velocity potential we must then have 

ip(r, 6)= A \n(r/a) + B6. (39) 

Note that, as we have previously hinted, the velocity potential is not a true 
function (because it has a discontinuity on going through 2n radians). The 
velocity potential must be interpreted as being defined patch- wise on overlapping 
regions surrounding the vortex core at r = 0. The velocity of the fluid flow is 

a^ Af + B °). (40) 

r 

Dropping a position-independent prefactor, the acoustic metric for a draining 
bathtub is explicitly given by 

ds 2 = -c 2 dt 2 + (dr - ^d?j +(rd9-^dt\ . (41) 

Equivalently 

ds 2 = - ( c 2 - ^ \ B2 \ dt 2 - 2- dr dt - 2B d9 dt + dr 2 + r 2 d9 2 . (42) 
\ r 2 J r 

A similiar metric, restricted to A = (no radial flow), and generalized to an 
anisotropic speed of sound, has been exhibited by Volovik |32| , that metric being 
a model for the acoustic geometry surrounding physical vortices in supcrfluid 
3 He. (For a survey of the many analogies and similarities between the physics 
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of superfluid 3He and the Standard Electroweak Model see Q , this reference is 
also useful as background to understanding the Lorentzian geometric aspects of 
3 He fluid flow.) Note that the metric given above is not identical to the metric 
of a spinning cosmic string, which would instead take the form |H] 

ds 2 = -c 2 (dt - Add) 2 + dr 2 + (1 - B)r 2 d6 2 . (43) 

In conformity with previous comments, the vortex fluid flow is seen to possess an 
acoustic metric that is stably causal and which does not involve closed timclikc 
curves. (At large distances it is possible to approximate the vortex geometry 
by a spinning cosmic string p^l , but this approximation becomes progressively 
worse as the core is approached.) 
The ergosphere forms at 



_ V^ 2 + 

i ergosphere ■ V 

Note that the sign of A is irrelevant in defining the ergosphere and ergo-region: 
it does not matter if the vortex core is a source or a sink. 

The acoustic event horizon forms once the radial component of the fluid 
velocity exceeds the speed of sound, that is at 

^horizon — • (^<-0 
C 

The sign of A now makes a difference. For A < we are dealing with a future 
acoustic horizon (acoustic black hole), while for A > we are dealing with a 
past event horizon (acoustic white hole). 

Though this construction has been phrased in (2 + 1) dimensions we are 
of course free to add an extra dimension by going to (3 + 1) dimensions and 
interpreting the result as a superposition of an ordinary vortex filament and a 
line source (or line sink). 

ds 2 = -c 2 dt 2 + (dr - -dt\ +frd9- —dt\ + dz 2 . (46) 



6 Slab geometries 

A popular model for the investigation of event horizons in the acoustic analogy 
is the one-dimensional slab geometry where the velocity is always along the z 
direction and the velocity profile depends only on z. The continuity equation 
then implies that p(z)v(z) is a constant, and the acoustic metric becomes 



ds 2 tx 



1 



v(z)c(z) 



-c{z) 2 dt 2 + {dz - v(z)dt} + dx 2 + dy 2 



(47) 
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That is 



ds 2 oc . , , , I- \c(z) 2 - v(z) 2 } dt 2 - 2v(z)dzdt + dx 2 + dy 2 + dz 2 ] . (48) 
v{zjc(z) 

If we set c = 1 and ignore the conformal factor we have the toy model acoustic 
geometry discussed by Unruh || page 2828, equation (8)] Jacobson ||, page 
7085, equation (4)], Corley and Jacobson ||, and Corley Q. (Since the con- 
formal factor is regular at the event horizon, we know that the surface gravity 
and Hawking temperature are independent of this conformal factor In the 

general case it is important to realise that the flow can go supersonic for either 
of two reasons: the fluid could speed up, or the speed of sound could decrease. 
When it comes to calculating the "surface gravity" both of these effects will 
have to be taken into account. 



7 The Painleve-Gullstrand line element 

To see how close the acoustic metric can get to reproducing the Schwarzschild 
geometry it is first useful to introduce one of the more exotic representations 
of the Schwarzschild geometry: the Painleve-Gullstrand line element, which is 
simply an unusual choice of coordinates on the Schwarzschild spacetime|] In 
modern notation the Schwarzschild geometry in ingoing (+) and outgoing (— ) 
Painleve-Gullstrand coordinates may be written as: 



ds 2 = -dt 2 + I dr ± J^^—dt ] + r 2 (d6 2 + sin 2 9 d<p 2 ) . (49) 



Equivalently 



ds 2 = - (l - dt 2 ± dr dt + dr 2 + r 2 (d6 2 + sin 2 9 dcf) . (50) 

This representation of the Schwarzschild geometry is not particularly well- 
known and has been rediscovered several times this century. See for instance 
Painleve jt5) , Gullstrand jl6) , Lemaitre jl7| , the related discussion by Israel Jl8| , 
and more recently, the paper by Kraus and Wilczek |ll| . The Painleve-Gullstrand 
coordinates are related to the more usual Schwarzschild coordinates by 



tpc — ts ± 
Or equivalently 



4M arctanh | J— — ) - 2 V 2GM r 



(51) 



j2GM/r , . 

dtpG = dts± T^GWr dr - (52) 



9 The Painleve-Gullstrand line element is often called the Lemaitre line element. 
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With these explicit forms in hand, it becomes an easy exercise to check the 
equivalence between the Painleve-Gullstrand line element and the more usual 
Schwarzschild form of the line clement. It should be noted that the + sign corre- 
sponds to a coordinate patch that covers the usual asymptotic region plus the re- 
gion containing the future singularity of the maximally extended Schwarzschild 
spacetime. It thus covers the future horizon and the black hole singularity. On 
the other hand the — sign corresponds to a coordinate patch that covers the 
usual asymptotic region plus the region containing the past singularity. It thus 
covers the past horizon and the white hole singularity. 

As emphasized by Kraus and Wilczek, the Painleve-Gullstrand line element 
exhibits a number of features of pedagogical interest. In particular the constant 
time spatial slices are completely flat — the curvature of space is zero, and all 
the spacetime curvature of the Schwarzschild geometry has been pushed into 
the time-time and time-space components of the metric. 

Given the Painleve-Gullstrand line clement, it might seem trivial to force 
the acoustic metric into this form: simply take p and c to be constants, and set 
v = i/2GM/r? While this certainly forces the acoustic metric into the Painleve- 
Gullstrand form the problem with this is that this assignment is incompatible 
with the continuity equation V • (pv) ^ that was used in deriving the acoustic 
equations. 

The best we can actually do is this: Pick the speed of sound c to be a 
position independent constant, which we normalize to unity (c = 1). Now 
set v = ^2GM/r, and use the continuity equation V ■ (pv) — to deduce 
p\v\ oc 1/r 2 so that p cx r~ 3 / 2 . Since the speed of sound is taken to be constant 
we can integrate the relation c 2 = dp/ dp to deduce the equation of state must 
be p = poo + c 2 p and that the background pressure satisfies p — Poo oc c 2 r~ 3 / 2 . 
Overall the acoustic metric is now 



ds 2 oc r~ 3/2 

The net result is conformal to the Painleve-Gullstrand form of the Schwarzschild 
geometry but not identical to it. For many purposes this is quite good enough: 
we have an event horizon, we can define surface gravity, we can analyze Hawking 
radiation. Since surface gravity and Hawking temperature are conformal invari- 
ants HU this is sufficient for analyzing basic features of the Hawking radiation 
process. The only way in which the conformal factor can influence the Hawk- 
ing radiation is through backscattering off the acoustic metric. (The phonons 
are minimally coupled scalars, not conformally coupled scalars so there will in 
general be effects on the frequency-dependent greybody factors.) 

If we focus attention on the region near the event horizon, the conformal 
factor can simply be taken to be a constant, and we can ignore all these com- 
plications. 
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—dt 



dr± 



2GM 



dt 



d9 2 



(53) 



8 The canonical acoustic black hole 



We can turn this argument around and ask: Given a spherically symmetric 
flow of incompressible fluid, what is the acoustic metric? What is the corre- 
sponding line element in Schwarzschild coordinates? If we start by assuming 
incompressibility and spherical symmetry, then since p is position independent 
the continuity equation implies v oc 1/r 2 . But if p is position independent then 
(because of the barotropic assumption) so is the pressure, and hence the speed 
of sound as well. So we can define a normalization constant ro and set 



v = c\ (54) 
The acoustic metric is therefore, up to an irrelevant position-independent factor, 

ds 2 = ~c 2 dt 2 + (dr ± c 4 dt) +r 2 {d0 2 + sin 2 9 d<j) 2 ). (55) 



If we make the coordinate change 

2/2 

dr = dt± M T ° \ . ... dr, (56) 

then 

dr 2 

ds 2 = -c 2 [l - {4/r 4 )}dr 2 + 4 + r 2 (d9 2 + sin 2 9 d<j) 2 ). (57) 

-1 { r o/ r ) 

This is not any of the standard geometries typically considered in general rela- 
tivity but is, in the sense described above, the canonical acoustic black hole. 

It is very important to realise that a time-dependent version of this canon- 
ical acoustic metric is very easy to set up experimentally p5[, since the time- 
dependent version of this canonical black hole metric is exactly the acoustic 
metric that is set up around a spherically-symmetric bubble with oscillating 
radius. For a bubble of radius R we have 




(58) 

We should only use this canonical metric for the fluid region outside the bubble, 
and only in the approximation that the ambient fluid is incompressible (e.g. 
water). For the typically gaseous and necessarily compressible medium inside 
the bubble (e.g. air) we should use a separate acoustic metric. The two acoustic 
metrics need not be continuous across the bubble wall. 

It is experimentally easy to generate (non-stationary) acoustic apparent hori- 
zons in this manner: In cavitating bubbles (typically air bubbles in water) it 
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is experimentally easy to get the bubble wall moving at supersonic speeds (up 
to Mach 10 in extreme cases). Once the bubble wall is moving supersonically 
an acoustic apparent horizon forms. It first forms at the bubble wall itself but 
then will typically detach itself from the bubble wall (since the apparent hori- 
zon will continue to be the surface at which the fluid achieves Mach 1) as the 
bubble wall goes supersonic. Since the bubble must eventually stop its collapse 
and re-expand, there is strictly speaking no acoustic event horizon (no absolute 
horizon) in this experimental situation, merely a temporary apparent horizon. 
(The apparent horizon must by construction last less than one sound-crossing 
time for the collapsing bubble.) 

To set up a geometry of this particular type with a true event horizon (or 
at the very least, an apparent horizon that lasts for many sound crossing times) 
requires a rather different physical setup: a big tank of fluid with a long thin 
pipe leading to the center. Then apply pressure to the tank till the outflow 
of fluid escaping through the pipe goes supersonic, being careful to maintain 
laminar flow and avoid turbulence. This would appear to be a technologically 
challenging project. 

9 Hawking radiation and "surface gravity" 

Establishing the existence of acoustic Hawking radiation follows directly from 
the original Hawking argument |36], [5?J once one realizes that the acoustic fluc- 
tuations effectively couple to the Lorentzian acoustic metric introduced above. 
The only subtlety arises in correctly identifying the "surface gravity" of an 
acoustic black hole. Because of the definition of event horizon in terms of 
phonons (null geodesies) that cannot escape the acoustic black hole, the event 
horizon is automatically a null surface, and the generators of the event horizon 
are automatically null geodesies. 

In the case of acoustics there is one particular parameterization of these 
null geodesies that is "most natural" , which is the parameterization in terms of 
the Newtonian time coordinate of the underlying physical metric. This allows 
us to unambiguously define a "surface gravity" even for non-stationary (time- 
dependent) acoustic event horizons, by calculating the extent to which this 
natural time parameter fails to be an affine parameter for the null generators 
of the horizon. (This part of the construction fails in general relativity where 
there is no universal natural time-coordinate unless there is a timelike Killing 
vector — this is why extending the notion of surface gravity to non-stationary 
geometries in general relativity is so difficult.) 

When it comes to explicitly calculating the surface gravity in terms of suit- 
able gradients of the fluid flow, it is nevertheless very useful to limit attention 
to situations of steady flow (so that the acoustic metric is stationary). This 
has the added bonus that for stationary geometries the notion of "acoustic sur- 
face gravity" in acoustics is unambiguously equivalent to the general relativity 
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definition. 

It is also useful to take cognizance of the fact that the situation simplifies 
considerably for static (as opposed to merely stationary) acoustic metrics. 



9.1 Static acoustic geometries 

To set up the appropriate framework, write the general stationary acoustic 
metric in the form 

ds 2 = - c [~c 2 dt 2 + (dx-v dt) 2 } . (59) 
The time translation Killing vector is simply — (1;0), with 

K 2 = 9liV K»K v ee -H^H 2 = - P -y - v 2 ]. (60) 

The metric can also be written as 

ds 2 = t \-{c 2 - v 2 )dt 2 -2v-dxdt + (dx) 2 ] . (61) 
c 

Now suppose that the vector v/(c 2 — v 2 ) is integrable, then we can define a new 
time coordinate by 

v • dx , 

dr = dt+^ =•. (62 

c — v z 

Substituting this back into the acoustic line element gives0 

,2 _ P 



ds' 

c 



(63) 



In this coordinate system the absence of the time-space cross-terms makes man- 
ifest that the acoustic geometry is in fact static (the Killing vector is hypersur- 
face orthogonal). The condition that an acoustic geometry be static, rather 
than merely stationary, is thus seen to be 



VX 1(^)1=°' ^ 

that is 

vxV(c 2 -v 2 ) = 0. (65) 

This requires the fluid flow to be parallel to another vector that is not quite the 
acceleration but is closely related to it. (Note that, because of the vorticity free 
assumption, |Vw 2 is just the three-acceleration of the fluid, it is the occurrence 
of a possibly position dependent speed of sound that complicates the above.) 

10 The corresponding formula in O is missing a factor of c and a bracket. 
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Once we have a static geometry, we can of course directly apply all of the 
standard tricks |3£| for calculating the surface gravity developed in general rela- 
tivity. We set up a system of fiducial observers (FIDOS) by properly normalizing 
the time-translation Killing vector 

Vfido 55 Fil = TWW^} (66) 

The four-acceleration of the FIDOS is defined as Afido = {Vfido • V)Vfido, 
and using the fact that if is a Killing vector, it may be computed in the standard 
manner 

AfIDO ~ + 2~\\K\F' (67) 

That is 

i r V ( c 2 -« 2 ) , v( P /c)l (6g) 



Afido — x 



(c 2 — v 2 ) {p/c) 

The surface gravity is now defined by taking the norm ||Afi.do||) multiplying 
by the lapse function, \\K\\ — y/ (p/c)\c 2 — v 2 }, and taking the limit as one 
approaches the horizon: |i>| — » c — remember this is the static case. The net 
result is 

UfidoW \\K\\ = ^v-V(c 2 -v 2 ) + 0(c 2 -v 2 ), (69) 
so that the surface gravity is given in terms of a normal derivative byp] 

1 d(c 2 - v 2 ) d(c-v) 

2 On on 

This is not quite Unruh's result 0, || since he implicitly took the speed of 
sound to be a position-independent constant. The fact that p drops out of 
the final result for the surface gravity can be justified by appeal to the known 
conformal invariance of the surface gravity [B4j]. Though derived in a totally 
different manner, this result is also compatible with the expression for "surface- 
gravity" obtained in the solid-state black holes of Reznik [Q , wherein a position 
dependent (and singular) refractive index plays a role analogous to the acoustic 
metric. As a further consistency check, one can go to the spherically symmetric 
case and check that this reproduces the results for "dirty black holes" enunciated 
in@. 

Since this is a static geometry, the relationship between the Hawking tem- 
perature and surface gravity may be verified in the usual fast-track manner — 
using the Wick rotation trick to analytically continue to Euclidean space p9| . 
If you don't like Euclidean signature techniques (which are in any case only ap- 
plicable to equilibrium situations) you should go back to the original Hawking 
derivations |3q, 13711 . 



11 Because of the background Minkowski metric there can be no possible confusion as to the 
definition of this normal derivative. 
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One final comment to wrap up this section: the coordinate transform we 
used to put the acoustic metric into the explicitly static form is perfectly good 
mathematics, and from the general relativity point of view is even a simplifica- 
tion. However, from the point of view of the underlying Newtonian physics of 
the fluid, this is a rather bizarre way of deliberately de-synchronizing your clocks 
to take a perfectly reasonable region — the boundary of the region of supersonic 
flow — and push it out to "time" plus infinity. From the fluid dynamics point 
of view this coordinate transformation is correct but perverse, and it is easier to 
keep a good grasp on the physics by staying with the original Newtonian time 
coordinate. 

9.2 Stationary but non-static acoustic geometries 

If the fluid flow does not satisfy the integrability condition which allows us 
to introduce an explicitly static coordinate system, then defining the surface 
gravity is a little trickier. The situation is somewhat worse than for general 
relativity since in the acoustic case we have no reason to believe that anything 
like the zeroth law of black hole mechanics holds [[lo) , nor do we have any reason 
to believe that stationary event horizons have to be Killing horizons. 

Recall that the zeroth law of black hole mechanics (constancy of the surface 
gravity over the event horizon) is proved in general relativity by appealing to 
the Einstein equations and imposing suitable energy conditions. In the acous- 
tic paradigm we have no analogue for the Einstein equations and no particular 
reason to suspect the existence of anything like a zeroth law. Sufficiently con- 
voluted supersonic flows would seem to be able to set up almost any pattern of 
surface gravity one wants. 

Similarly, in general relativity the fact that stationary but non-static black 
holes possess Killing horizons is related to the axisymmetry that is deduced 
from the fact that non-axisymmetric black holes are expected to lose energy via 
gravitational radiation and so dynamically relax to an axisymmetric configura- 
tion — in the fluid dynamic models discussed here I have explicitly allowed for 
external driving forces and explicitly excluded back reaction effects, therefore 
there is no particular reason to expect acoustic black holes to dynamically re- 
lax to axisymmetry. In particular, this means that even for stationary acoustic 
geometries there is no particular reason to expect the acoustic event horizon to 
in general be a Killing horizon. 

So what does survive of our usual general relativistic notions for acoustic 
event horizons in stationary but non-static geometries? Recall that by con- 
struction the acoustic apparent horizon is in general defined to be a two-surface 
for which the normal component of the fluid velocity is everywhere equal to the 
local speed of sound, whereas the acoustic event horizon is characterized by the 
boundary of those null geodesies (phonons) that do not escape to infinity. In 
the stationary case these notions coincide, and it is still true that the horizon is 
a null surface, and that the horizon can be ruled by an appropriate set of null 
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curves. Suppose we have somehow isolated the location of the acoustic horizon, 
then in the vicinity of the horizon we can split up the fluid flow into normal and 
tangential components 



Vj_ 



where 



(71) 



Here (and for the rest of this particular section) it is essential that we use the 
natural Newtonian time coordinate inherited from the background Newtonian 
physics of the fluid. In addition n is a unit vector field that at the horizon is 
perpendicular to it, and away from the horizon is some suitable smooth exten- 
sion. (For example, take the geodesic distance to the horizon and consider its 
gradient.) We only need this decomposition to hold in some open set encom- 
passing the horizon and do not need to have a global decomposition of this type 
available. Furthermore, by definition we know that v± = c at the horizon. Now 
consider the vector field 

L" = (l;w|). (72) 

Since the spatial components of this vector field are by definition tangent to a 
constant time slice through the horizon, the integral curves of this vector field 
will be generators for the horizon. Furthermore the norm of this vector (in the 
acoustic metric) is 



\L\? = -^ 



— (c 2 — v 2 ) — 2uii • v + vu ■ «|| cx (c 2 — Vj_) 



(73) 



In particular, on the acoustic horizon L M defines a null vector field, the integral 
curves of which are generators for the acoustic horizon. I shall now verify that 
these generators are geodesies, though the vector field L is not normalized with 
an affine parameter, and in this way shall calculate the surface gravity. (For 
clarity, I will drop the conformal factor because I already know that it will not 
affect the surface gravity pif.) 

Consider the quantity (L ■ V)L and calculate 



L a V a L» = L a (V a L p 

To calculate the first term note that 

P i 



V L a )g^+±Vp(L 2 )g^. 



(74) 



(75) 



Thus 



K/3] 



And 



L a L 



1PM 



-v, [£(< 



-V, [f(c 2 



c 



vi)) 



hJl 



(76) 



(77) 
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On the horizon, where c = v±, this simplifies tremendously 

(L a L [0M )\ horlzon = -t (OjV^c 2 -vl)). (78) 
Similarly, for the second term we have 

Vp{L*)=(Q-V J \Py-vl)\). (79) 



On the horizon this again simplifies 

V p(L 2 )\ ho „ zon = +^ (0;V,(c 2 -vl)). (80) 
There is partial cancellation between the two terms, and so 

L ° VaL " = + h (^iKc 2 - «i)]5 - uV'JViKc 3 - vl)}) . (81) 

But, as we have already seen, at the horizon the gradient term is purely normal. 
Thus 

L^ = + _-L gr ii(l; W {). (82) 
Comparing this with the standard definition of surface gravity p(|P| 



we finally have 



L a W a L" = +^-L^. (83) 

l d(c 2 -vp d(c-v ± ) 
9H =2 dn = °—d^- (84) 
This is in agreement with the previous calculation for static acoustic black holes, 
and insofar as there is overlap, is also consistent with results of Unruh ||, ||, 
Reznik [jl4j , and the results for "dirty black holes" (2TJ] . From the construction 
it is clear that the surface gravity is a measure of the extent to which the 
Newtonian time parameter inherited from the underlying fluid dynamics fails 
to be an affine parameter for the null geodesies on the horizon. 



10 Geometric Acoustics 

Up to now, have been developing general machinery to force acoustics into 
Lorentzian form. This can be justified either with a view to using fluid me- 
chanics to teach us more about general relativity, or to using the techniques of 
Lorentzian geometry to teach us more about fluid mechanics. 

12 There is an issue of normalization here. On the one hand we want to be as close as 
possible to general relativistic conventions. On the other hand, we would like the surface 
gravity to really have the dimensions of an acceleration. The convention adopted here is the 
best compromise I have come up with. 
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For example, given the machinery developed so far, taking the short wave- 
length/high frequency limit to obtain geometrical acoustics is now easy. Sound 
rays (phonons) follow the null geodesies of the acoustic metric. Compare this to 
general relativity where in the geometrical optics approximation light rays (pho- 
tons) follow null geodesies of the physical spacetime metric. Since null geodesies 
are insensitive to any overall conformal factor in the metric [^8[ ^9], [30| one might 
as well simplify life by considering a modified conformally related metric 



(85) 



This immediately implies that, in the geometric acoustics limit, sound propaga- 
tion is insensitive to the density of the fluid. In this limit, acoustic propagation 
depends only on the local speed of sound and the velocity of the fluid. It is only 
for specifically wave related properties that the density of the medium becomes 
important. 

We can rephrase this in a language more familiar to the acoustics community 
by invoking the Eikonal approximation. Express the linearized velocity poten- 
tial, in terms of an amplitude, a, and phase, ip, by tp\ ~ ae ltp . Then, ne- 
glecting variations in the amplitude a, the wave equation reduces to the Eikonal 
equation 

hf" d^ip d v ip = 0. (86) 

This Eikonal equation is blatantly insensitive to any overall multiplicative pref- 
actor (conformal factor). 

As a sanity check on the formalism, it is useful to re-derive some standard 
results. For example, let the null geodesic be parameterized by x^it) = (t; x(t)). 
Then the null condition implies 



dx 11 dx v 




— r- = 


dt 


dt 








dx 




-dt~ V0 



dx 1 
] ~dt 



dx 1 dx 1 
~~dt~dt 



= 



(87) 



Here the norm is taken in the flat physical metric. This has the obvious inter- 
pretation that the ray travels at the speed of sound, c, relative to the moving 
medium. 

Furthermore, if the geometry is stationary one can do slightly better. Let 
x^(s) = (t(s); x(s)) be some null path from x\ to £2, parameterized in terms of 
physical arc length (i.e. ||<&r/<fs|| = 1). Then the tangent vector to the path is 



dxf* _ (dt dx 1 
ds yds 1 ds 



(88) 
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The condition for the path to be null (though not yet necessarily a null geodesic) 
is 

dx^ dx u n 

^ir* =0 - (89) 

Using the explicit algebraic form for the metric, this can be expanded to show 



dt \ 2 _ . ( dx 1 \ ( dt 



Solving this quadratic 



) = V 1 \ -2 -• (91) 



ds) 

Therefore, the total time taken to traverse the path is 



T[i\ = I (dt/ds) ds 

-i— j{ J(c 2 - vl)ds 2 + (vj dx 1 ) 2 - vl dx 1 ]. (92) 



c 2 

If we now recall that extremizing the total time taken is Fermat's principle for 
sound rays, we see that we have checked the formalism for stationary geome- 
tries (steady flow) by reproducing the discussion on page 262 of Landau and 
Lifshitz §|. 

As a second example of the insights arising from the Lorentzian point of view 
consider the "reciprocity theorem" . Suppose a pulse of sound is emitted at time 
t\ at position x\. The disturbance propagates according to the inhomogeneous 
differential equation 



AV> = -^=S 4 (x^ - a£) = -^5{t- h) 5 3 (x - xx). (93) 



-9 P 
The solution to this is the retarded scalar Green function 

source at x\ 

= G R (x,xi). (94) 

The Green function has well known symmetry properties that are completely 
unaffected by any time dependence in the underlying acoustic metric. We may 
in the usual manner, construct advanced and retarded Green functions that 
vanish outside the past and future sound cones respectively. Then 

G R (x 2l xi) = Ga(xi,x 2 ). (95) 

So that the reciprocity theorem for the velocity potential is valid in absolute 
generality. 

^X2)\source at „ = ) I foZrclZlf ■ (96) 
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To get a reciprocity theorem for pressure one has to recall 



Pi = Po{dti>\ + v ■ VV>i). 



(97) 



Then, by restricting to the case of fluid at rest (vq = 0, dtpo = 0, dtpo = 0), using 
the time translation invariance of the Green functions, and the time reversal 
property Ga Gr, one has 



source at x\ 



Pi 
Pa 



(#1,^2 ^l) | source at 



(98) 



This result is still much more general than the usual reciprocity theorem. 



11 Limitations 

The derivation of the wave equation made two key assumptions — that the flow 
is irrotational flow and the fluid is barotropic. 

The d'Alembertian equation of motion for acoustic disturbances, though 
derived only under the assumption of irrotational flow,^| continues to make 
perfectly good sense in its own right if the background velocity field vq is given 
some vorticity. This leads one to hope that it might be possible to find a suitable 
generalization of the present derivation that might work for flows with nonzero 
vorticity. In this regard, note that if the vorticity is everywhere confined to thin 
vortex filaments, the present derivation already works everywhere outside the 
vortex filaments themselves. 

The technical problem with flows with non-zero vorticity is that the vorticity 
in the background flow couples to the perturbations and generates vorticity 
in the fluctuations. Then sound waves can no longer be represented simply 
by a scalar potential and a much more complicated mathematical structure 
results. (Phonons are no longer simply minimally coupled scalar fields and the 
appropriate generalization is sufficiently unpleasant as to be intractable.) 

The restriction to a barotropic fluid (p a function of p only) is in fact also 
related to issues of vorticity. Examples of barotropic fluids are: 

• Isothermal fluids subject to isothermal perturbations. 

• Fluids in convective equilibrium subject to adiabatic perturbations. 

See for example §l|, §311, pp 547-548, and §313 pp 554-556. Failure of the 
barotropic condition implies that the perturbations cannot be vorticity free and 
thus requires more sophisticated analysis. 

If the fluid is in addition inviscid then the analysis of this paper implies 
a hidden Lorentz invariance in the acoustic equations. This hidden Lorentz 

13 But remember that irrotational flow is automatic for superfluids [ ^t , and is natural in 
situations of high symmetry. 
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invariance is more than just a formal quirk: If one has "atoms" held together 
by phonons (Cooper pairs?), then these atoms, and complex systems built up 
out of such atoms, will see (hear) an acoustic special relativity that is as real 
to them as Einstein's special relativity is to us. Furthermore these systems 
would with additional observation detect (hear) an acoustic general relativity 
- but instead of the Einstein-Hilbert equations of our general relativity they 
would experience an acoustic general relativity governed by the hydrodynamic 
equations. 

If the fluid has non-zero viscosity then there will be violations of this acoustic 
Lorentz symmetry. These violations are momentum dependent and, as I shall 
discuss in the next section, they are small at low momentum. 



12 Viscosity: breaking the Lorentz symmetry 

After this long build-up emphasizing the hidden Lorentzian geometry hiding in 
(inviscid vorticity-free barotropic) fluid dynamical equations, I will now show 
how to explicitly break the Lorentz symmetry. From the atomic perspective 
underlying continuum fluid mechanics the eventual breakdown of the Lorentz 
symmetry governing the notion of the phonons is no great surprise: eventually, 
once the wavelength of the phonons is less than the mean interatomic spacing 
in the fluid, we should certainly expect modifications to the phonon dispersion 
relation jl], [j). Specific ad hoc mutilations of the dispersion relation have been 
considered by Jacobson |l| , Unruh Q , Corley and Jacobson |)| , and Corley |ll| . 
I shall now show that a similar but not identical breakdown of acoustic Lorentz 
invariance can be deduced directly from the continuum equations merely by 
adding the effects of viscosity. 

Of course, the fundamental equations of fluid dynamics, the equation of 
continuity @ and Euler's equation @ are unaltered. What changes is the 
expression for the driving force in Euler's equation so that (Q) becomes [ ^H 
§328, page 576-577] 

F = -Vp - pV(/> - + pi> (\7 2 v + ^V(V • v)j . (99) 

Here v denotes kinematic viscosity. I again take the flow to be vorticity free, 
and again choose the fluid to be barotropic. Repeating the steps that led to ( |To| ) 
now show that Euler's equation reduces to 

- d t ij) + h + i(V?/>) 2 + (p + $ + -v V 2 ip = 0. (100) 
£ o 

This again is a well-known equation, simply being Burgers' equation subject to 
external driving forces |4l|] . (In obtaining this equation it is necessary to assume 
that the kinematic viscosity v is position independent. In addition it is common 
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Rearranging 



practice, though not universal, to absorb the 4/3 into a modified definition of 
kinematic viscosity.) 

Linearization proceeds as previously. For the continuity equation there are 
no changes, while linearizing the Euler equation (Burgers' equation) yields 

-dtih + h + i(W> ) 2 + + $ + ^ vVo = 0. (101) 

-dtipi + — -v -VTpi + -v VVi = 0. (102) 
Po 3 

Pi = Po (d t ipi + v ■ Vipi - -fV 2 ^ ■ (103) 

As before, we substitute this linearized Euler equation into the linearized con- 
tinuity equation, to obtain the physical wave equation: 

- d t po (dtip! + v ■ Vip! - |i/ V 2 ^] 

Po V^i - ^ po v (dtip! + vo ■ W>i - V 2 Vi^ = 0.(104) 

Using the same matrix f^ u defined previously the above wave equation is easily 
rewritten as[j^| 

W" Wi) = ~\po v (jfi + v ■ V\ [c- 2 V 2 ^] . (105) 
In terms of the d'Alembertian associated with the acoustic metric this reads 

A^ = -t-(|- + €b ■ v) [c- 2 V 2 ^] . (106) 



3 po \9t 

The convective derivative appearing here may easily be converted into four- 
dimensional form by utilizing the acoustic four-velocity for the fluid. Recall 
that 

V» = (107) 
y/PoC 

It is easy to see that this is a timelike unit vector in the acoustic metric, so that 

4 uc^ 

A^i - -r } ^= (V%) [c- 2 V 2 ^] . (108) 
3 ^fpoc 

The V 2 ^i term explicitly couples only to the flat spatial metric and can be 
written in terms of the acoustic metric by noting that 

<T = -VT + - ^g^ace- (109) 



14 I have used the continuity equation for the background fluid flow to pull the factor po 
outside the convective derivative. 
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It is the explicit appearance of the fluid four-velocity in the above expressions 
that justifies my claim that viscosity breaks the acoustic Lorentz invariance. 



Sanity check I: If the background fluid flow is at rest and homogeneous 
(vq = 0, and with po and c independent of position) then this viscous wave 
equation reduces to 

dfVi =c 2 V 2 Vi + ^5 t V 2 Vi. (HO) 
This equation may be found, for instance, in §359 pages 646-648 of Lamb pl| . 

Sanity check II: Take the Eikonal approximation in the form 

tpi = a(x) exp(— i[u>t — k • x\ ), (m) 

with a(x) a slowly varying function of position. Furthermore, agree to ignore 
derivatives of the metric. Then the viscous wave equation in the Eikonal ap- 
proximation reduces to 

- (u-v- kf + c 2 k 2 - iv-{w-v ■ k)k 2 = 0. (112) 

o 

This lets us write down a dispersion relation for sound waves 



- r o,o f2vk 2 \ 2 2vk 2 , . 

u = v ■ k± Jc 2 k 2 - f — — j (U3) 

The first term simply arises from the bulk motion of the fluid. The second 
term specifically introduces dispersion due to viscosity, while the third term is 
specifically dissipative. The ad hoc models introduced in Jacobson j3j, Unruh B, 
Corley and Jacobson |1, and Corley |ll| are exactly recovered by ignoring the 
dissipation due to viscosity but retaining the dispersion due to viscosity. 

Note that the violation of Lorentz invariance is suppressed at low momentum. 
This is in agreement with general arguments of Nielsen et al jl2| flij , though 
it should be borne in mind that Nielsen et al were dealing with interacting 
quantum field theories and the context here is, if not purely classical, at worst 
one of free phonons propagating on a fixed classical background. (An alternative 
model for the breakdown of Lorentz invariance has been discussed by Everett [^5| , 
(46f.) The violations of Lorentz symmetry become significant once 

fc«fc = -. (114) 
v 

But from the atomic theory of (normal) fluids 
mean free path 



mean free time 



c • mean free path. (H5) 
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This gives the very reasonable result 

mean free path' ^ 

verifying that the breakdown of acoustic Lorentz invariance is explicitly linked 
to the atomic nature of matter. 



13 Precursors 



It is perhaps surprising that anything new can be said about so venerable a 
subject as fluid dynamics. Certainly there are precursors to the discussion of 
this paper in the fluid dynamics literature. For instance, take the background 
to be static, so that vq = 0, while dtpo = = dtpo, though po and hence c 
are permitted to retain arbitrary spatial dependencies. Then the wave equation 
derived in this paper reduces to 

9 t V = c 2 -V • (p V# (H7) 
Po 

This equation is in fact well known. It is equivalent, for instance to eq. (13) 
of §313 of Lamb's classic Hydrodynamics f2l|| . See also eq. (1.4.5) of the recent 
book by DeSanto |}47|. The superficially similar wave equations discussed by 
Landau and Lifshitz j§§ (see §74, eq. (74.1)), and by Skudrzyk || (see p 282), 
utilize somewhat different physical assumptions concerning the behaviour of the 
fluid. 

In a somewhat different vein, the modern study of classical continuum me- 
chanics has greatly benefited from the use of 3-dimensional Riemannian geom- 
etry to describe the physics of the spatial configurations of elastic media and 
other continua |l8| ^FJ, |5(J. Analyses of this type have traditionally treated 
space and time on quite separate footings. 

The most direct precursor of the results derived in this paper are due to 
Unruh [1 and Jacobso n Q , and in the body of work prompted by those pa- 
pers |, 5, 1 1, |, | M In, [l| Fl p. 



14 Summary and Discussion 

Acoustic waves in an inviscid fluid can, under the assumptions of irrotational 
barotropic flow, be described by an equation of motion involving the scalar 
d'Alembertian of a suitable Lorentzian geometry. For inhomogeneous flows this 
Lorentzian geometry will exhibit nonzero Riemann curvature. 

Traditionally, Lorentzian geometries have been of interest to physics only 
within the confines of Einstein's theory of gravitation. The results of this paper 
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provides the general relativity community with a very down to earth physi- 
cal model for certain classes of Lorentzian geometry. This is of interest both 
pedagogically and because it extends the usefulness of Lorentzian differential 
geometry beyond the confines of Einstein gravity. 

Particularly intriguing is the fact that while the underlying physics of fluid 
dynamics is completely nonrelativistic, Newtonian, and sharply separates the 
notions of space and time, one nevertheless sees that the acoustic fluctuations 
couple to a full-fledged Lorentzian spacetime. 

As discussed by Unruh , (and subsequent papers ||, ||, [|, H @, H H> 
E2I |l3|, an acoustic event horizon will emit Hawking radiation in the form 
of a thermal bath of phonons at a temperature 

kT H = (118) 

Z7TC 

(Yes, this really is the speed of sound in the above equation, and is re- 
ally normalized to have the dimensions of a physical acceleration.) Using the 
numerical expression 

T H = (1.2 x lcr 9 ^™) 

it is clear that experimental verification of this acoustic Hawking effect will 
be rather difficult. (Though, as Unruh has pointed out |l]], this is certainly 
technologically easier than building [general relativistic] micro-black holes in 
the laboratory.) 

A particularly important side effect of this entire analysis is that it forces us 
to re-examine all of black hole physics to cleanly separate what is intrinsic to 
general relativity from what is generic to Lorentzian geometries. The acoustic 
analog for black hole physics accurately reflects half of general relativity — the 
kinematics due to the fact that general relativity takes place in a Lorentzian 
spacetime. The aspect of general relativity that does not carry over to the 
acoustic model is the dynamics — the Einstein equations. Thus the acoustic 
model provides a very concrete and specific model for separating the kinematic 
aspects of general relativity from the dynamic aspects. 

In particular, perhaps the most important lesson to be learned is this: Hawk- 
ing radiation from event horizons is a purely kinematic effect that occurs in any 
Lorentzian geometry with an event horizon and is independent of any dynam- 
ical equations imposed on the Lorentzian geometry. On the other hand, the 
classical laws of black hole mechanics |Q are intrinsically results of the dynam- 
ical equations (Einstein equations) that have no analog in the acoustic model. 
Thus Hawking radiation persists even in the absence of the laws of black hole 
mechanics and, in particular, the existence or otherwise of Hawking radiation is 
now seen to be divorced from the issue of the existence or otherwise of the laws 
of black hole thermodynamics. Hawking radiation is a purely kinematical effect 



mnnms- 



1 o{c — v± ) 
r. Bn. 



(119) 
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that will be there regardless of whether or not it makes any sense to assign an 
entropy to the event horizon — and attempts at deriving black hole entropy 
from the Hawking radiation phenomenon are thereby seen to require specific 
dynamical assumptions about the (at least approximate) applicability of the 
Einstein equations. 
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